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Biases for neutron-star mass, radius and distance 
measurements from Eddington-limited X-ray bursts 
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ABSTRACT 

Eddington-limited X-ray bursts from neutron stars can be used in conjunction with 
other spectroscopic observations to measure neutron star masses, radii, and distances. 
In order to quantify some of the uncertainties in the determination of the Eddington 
limit, we analysed a large sample of photospheric radius-expansion thermonuclear 
bursts observed with the Rossi X-ray Timing Explorer. We identified the instant at 
which the expanded photosphere "touches down" back onto the surface of the neutron 
star and compared the corresponding touchdown flux to the peak flux of each burst. 
We found that for the majority of sources, the ratio of these fluxes is smaller than 
~ 1.6, which is the maximum value expected from the changing gravitational redshift 
during the radius expansion episodes (for a 2M Q neutron star). The only sources for 
which this ratio is larger than ~ 1.6 are high inclination sources that include dippers 
and Cyg X-2. We discuss two possible geometric interpretations of this effect and show 
that the inferred masses and radii of neutron stars are not affected by this bias. On 
the other hand, systematic uncertainties as large as ~ 50% may be introduced to the 
distance determination. 
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1 INTRODUCTION 

Photospheric radius-expansion (PRE) thermonuclear bursts 
are a key observational tool for studies of low-mass X- 
ray binaries. These systems consist of a neutron star ac- 
creting from a low-mass stellar companion. Thermonuclear 
(type I) bursts arise when a critical temperature and den- 
sity is reached in the accumulated fuel layer, and unsta- 
ble H or He ignition takes place (e.g. iLewin et al.1 Il993l ; 
IStrohmaver fc Bildstenl [2006). The accumulated fuel then 
burns and is exhausted within 10-100 s. Ongoing accretion 
leads to subsequent bursts separated by hours to tens of 
hours (depending upon the accretion rate). 

The X-ray flux early in the burst can be sufficient to ex- 
ceed the local Eddington limit at the surface of the neutron 
star. In that case, radiation presure drives an expansion of 
the photosphere, and excess flux is converted to kinetic and 
gravitational potential energy. The expansion of the pho- 
tosphere at approximately constant luminosity produces a 
characteristic pattern of the X-ray spectral variation in these 
bursts, giving rise to a peak in the blackbody radius at the 
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same time as a local minimum in the blackbody tempera- 
ture. Measurement of the X-ray flux during the PRE episode 
allows determination of the distance to the burst source (e.g. 
Basinska et al. 1984 ). as well as the mass and radius of the 

1 " 1 1 1| — , p-n 1| —j 

neutron star fe.g. lDamen et al . 1990; Ozel 200a) These mea- 
surements are hampered by intrinsic variation in the Ed- 
dington lumi nosities between sour ces and even from burst 
to bu rst (e.g. iKuulkers et afll2003l ; see also iGallowav et all 
I2008D . 

One subtlety that arises for measurements of the Ed- 
dington flux -pEdd is the precise instance during the burst at 
which this value is measured. The observed flux during the 
PRE episode is expected to vary due to the changing grav- 
itational redshift resulting from the varying height of the 
photosphere above the neutron star. Ideally, the flux should 
be measured when the Eddington limit is first reached, 
just before radius expansion commences, or alternatively 
when radius expansion ceases and the expanded photosphere 
"touches down" once again onto the neutron star. The in- 
stance when radius expansion commences is typically diffi- 
cult to identify, since it occurs during the burst rise. Instead, 
the flux at touchdown, at the end of the PRE episode, has 
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long been th ought to be the bes t time to measure the Ed- 
dington flux (|Damen et al.lll990h . 

Recently, measurements of the X-ray flux at touch- 
down in radius-expansion bursts from EXO 0748—676 were 
used (along with other observational constraints) to de- 
rive limits on the mass and radius of the neutron star 
ijOzell l2006h . The Eddington flux at touchdown measured 
from bursts observed by EXOSAT and RXTE, 2.25 x 
10 -8 ergcm _2 s _1 , is significantly lower than the maxi- 
mum flux reported fro m the bursts obse rved by RXTE, of 
5.2 x 10" 8 erg c m" 2 s" 1 dWolff et al]|2005h . This discrepancy 
(also noted by iDamen et al.l 1990t) of more than a factor of 
2, is far in excess of that expected from the varying gravi- 
tational redshift during radius expansion. If the touchdown 
flux does indeed correspond to the Eddington limit, it is dif- 
ficult to understand how the flux earlier in the burst could 
exceed that value by such a large margin. 

Motivated by this result, we undertook a broader study 
of the relationship between the peak radius-expansion burst 
flux and the flux at touchdown in different burst sources. In 
Sj2]we describe the observations and analysis of the bursts. 
In Sj3] we present the results of our study. In section £|3 . 1 1 
we identify sources with low touchdown fluxes (compared to 
the peak fluxes), while in i !3.2l we investigate the effect of the 
maximum radius reached during the expansion. We present 
the broader conclusions of this study in Sj4] 



minimum in fcTbb- By this criteria, 246 bursts in the cata- 
logue exhibited unambiguous evidence of radius expansion. 

We subsequently analysed the bolometric flux evolu- 
tion of each of these bursts to determine the touchdown 
flux Ft. Immediately prior to touchdown, the photosphere 
is expected to be contracting back towards the neutron 
star, reflected by the decreasing blackbody radius. The long 
(~ 1 s) observed timescale of the radius contraction (com- 
pared to the ~ 1 ms free-fall timescale) suggests that the 
atmosphere is still primarily radiation-pressure supported 
during the contraction phase. Because of this condition, the 
flux remains equal to the local Eddington limit, and the 
blackbody temperature must be increasing. Immediately af- 
ter touchdown, the photospheric radius is roughly constant, 
but now the burning has more or less ceased and the neutron 
star gradually cools, giving a decreasing flux and blackbody 
temperature. We thus identify the time at which the pho- 
tosphere touches down as the time of a local maximum in 
the blackbody temperature, following the time of maximum 
radius. 

This definition is illustrated in the two examples of Fig. 
[1] The time of touchdown can be unambiguously determined 
from the time-history of fcTbb; in both cases a clear local 
maximum can be seen (middle panels), following the radius 
maximum (bottom panels). While not all the cases were as 
clear-cut as these examples, we nevertheless were able to 
determine the time of touchdown for each of the radius- 
expansion bursts in the sample. 



2 OBSERVATIONS AND ANALYSIS 

We used analyses of thermonuclear bursts observed by the 
R ossi X-ray Timing Exp lorer (RXTE) from the catalogue 
of iGallowav et al. (2008). This catalogue contains analyses 
of all bursts present in public data since shortly after the 
RXTE launch on 1995 December 30. The primary instru- 
ment f or burst studies is the Proportional Counter Array 
fPCA; [jahoda et al.lll996T ). consisting of five identical, co- 
aligned proportional counter units (PCUs), sensitive to pho- 
tons in the energy range 2-60 keV and with a total effective 
area of ~ 6500 cm 2 . 

Data were analysed with LHEASOFT version 5.3 (re- 
leased 2003 November 17) or later; this release was the 
first to incorporate a significant reduction of the geomet- 
ric area of the PCUs for improved consistency with cal- 
ibration sources (| Jahoda et al.ll2006h . Data with the best 
spectral and temporal resolution were used to extract time- 
resolved spectra covering the PCA bandpass every 0.25 s 
near the peak of each burst, becoming longer in the burst 
tail to preserve signal-to-noise. Each burst spectrum was 
fitted using an absorbed blackbody, with a spectrum ex- 
tracted from a (typically) 16 s interval prior to the burst 
as the backgro und, following standard X-ray burst analysis 
proce dure (e.g. Ivan Paradiis fc Lewinl ll986; Ku ulkers et al.l 
l2002i ). To take into account gain variations over the life of 
the instrument, we generated a separate response matrix for 
each burst using PCARSP version 10.1. 

The bolometric flux at each timestep was calculated 
from the fitted blackbody temperature fcTbb and radius Rbb- 
Each burst was inspected visually for evidence of radius ex- 
pansion, specifically a local maximum in i?bb near the time 
the maximum burst flux was reached, synchronous with a 



3 RESULTS 

We measured the flux at touchdown Ft for each of 246 pho- 
tospheric radius-expansion bursts in our sample, and calcu- 
lated its ratio to the peak flux / = F p /F t . For 61 bursts 
(around 1/4 of the sample) the maximum flux was achieved 
at touchdown, so that F v — F t and / = 1. More gen- 
erally, F t was consistent with F p to within the 3<r confi- 
dence limit (based on the uncertainties of F p and Ft) for 
184 of the bursts, or almost 3/4 of our sample (the burst 
shown in the left panels of Fig. [T]is such an example). In 
many cases this commensurability of the peak and touch- 
down fluxes arises from a common behaviour in PRE bursts, 
where the flux continues to increase following the time the 
maximum radius is r eached (as noted for 4U 1636—536 by 
ISugimoto et al.1 [1984T ). We note that / > 1 for all three 
radius-expansion bursts from EXO 0748—676, as we e xpect 
based on t he comparison of the published fluxes by I6zell 
(|2006h and IWolff et al.l l|2005l ). Thus, we made a closer ex- 
amination of what conditions give rise to values of / > 1. 

3.1 Sources with low touchdown fluxes 

Radius-expansion bursts from eighteen different sources 
gave rise to values of / > 1. In many cases these values were 
only a few percent in excess of 1; by comparison, the /-value 
for two of the three bursts from EXO 0748—676 were in ex- 
cess of / > 1.9, and the maximum value over all the bursts 
was / = 6.72. Bursts with / > 1.6 arose from just seven 
sources: EXO 0748-676, MXB 1659-298, 4U 1916-053, 
GRS 1747-312, XTE J1710-281, 4U 1746-37 and Cyg X-2. 
The remarkable feature of this selection on / is that six of the 
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Figure 1. Example photospheric radius-expansion bursts from 
4U 1636-536 (left panels) and EXO 0748-676 (right panels). 
From top to bottom, the panels in each column show the bolo- 
metric luminosity (units of 10 38 erg s — 1 ); blackbody temper- 
ature (keV); and blackbody radius (km). The luminosity and 
blackbod y radii are calculate d for a distance to 4U 1636—536 
of 6 kpc [Galloway et alll2006r) , and to EXO 0748-676 of 7.4 kpc 
llGallowav et al.ll200Sir The open circles indicate the time of peak 
burst flux, while the open diamonds indicate the time when the 
photosphere is assumed to have just returned to the neutron star 
surface (the "touchdown" point). Note that the flux at touchdown 
Ft for 4U 1636—536 is essentially identical to the maximum flux, 
while for EXO 0748—676 it is approximately half that value. 



Figure 2. The distribution of / = F p /Ft, the ratio of the peak 
flux to the touchdown flux, for 246 radius-expansion bursts in 
our sample. The sample has been divided based on the source 
type; we plot the distribution for the dippers (EXO 0748—676, 
MXB 1659-298, 4U 1916-053, GRS 1747-312, XTE J1710-281 
and 4U 1746—37, solid histogram) separately from the non- 
dipping sources (dotted histogram). The overwhelming majority 
of bursts from the non-dipping sources have / < 1.2; we have 
indicated this in the histogram by labelling the first bin (which 
is cut off) with the data value. The inset shows the three bursts 
with / > 4 at the extreme end of the / distribution. Note that 
the only non-dipping source producing bursts with / > 1.6 is 
Cyg X-2. 



seven sources listed above (excluding Cyg X-2) all exhibit 
regular X-ray dips (see e.g. iLiu et all 1200 ll ) attributable to 
absorption from structures in the accretion disk, occuring at 
preferential orbital phases. Since the accretion disk is usu- 
ally thought to have limited vertical extent, the presence of 
dips is normally taken as an indication of high system incli- 
nation, i.e., i ~ 90°. We thus conclude that systematically 
large /-values overwhelmingly (but not exclusively) occur 
in bursts from systems with high inclination; we discuss the 
remaining example, Cyg X-2, in fJJ] 

We thus divided the sample between dipping and non- 
dipping sources, and calculated the combined distribution 
of / from each sample (Fig. [5]). While there is substantial 
overlap between the two distributions, the preponderance of 
bursts with / ~ 1 for the non-dipping sources makes the 
two distributions highly discrepant. The K-S test value for 
the two samples is 0.886, with associated probability 2 x 
10~ 21 ; this indicates unambiguously that the two samples 
arise from different populations. For the dippers, the mean 
/ = 2.1 ±0.9, so that the flux at touchdown was, on average, 
less than half the maximum flux. 



3.2 The effect of the maximum radius in dipping 
sources 

While the association of low touchdown fluxes with high 
system inclination is clear, the underlying cause (as well 
as the wide range of the / values for the dipping sources) 
begs an explanation. The large /-range was most notable for 
MXB 1659—298, where the /-values were between 1.3 and 
5.1 for the 12 radius-expansion bursts in the sample. 

The high end of the distribution in the ratio / is too 
large to be explained by the change in the gravitational red- 
shift during the expansion of the photosphere. For example, 
for a 10 km, 2Mq neutron star, the flux ratio arising from 
redshift changes can be at most equal to 1.6. Instead, the 
clear separation between the /-distributions of the dippers 
and non-dippers strongly suggests a geometric origin of the 
/-ratios. 

A geometric interpretation is further supported by the 
correlation between the ratio / observed in dippers and 
the maximum radius attained by the photosphere during 
the radius expansion episode (Figure [3} . In order to calcu- 
late the maximum radius, we corrected for the distance to 
the dipping sources based (in the case of EXO 0748—676, 
MXB 1659-298, 4U 1916-053 and XTE J1710-281) on 
the distance derived f rom the peak flux of the PRE bursts 
jGallowav et al.ll2008l ). In the case of GRS 1747-312 and 
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Figure 3. The ratio of peak to touchdown flux / = Fp/Ft 
as a function of the derived maximum photospheric expan- 
sion radius for a sample of bursts observed from 6 dip- 
ping sources: EXO 0748-676, MXB 1659-298, 4U 1916-053, 
XTE J1710-281, GRS 1747-312 and 4U 1746-37. The dashed 
lines shows the best linear fit excluding the single outlier, an ex- 
tremely intense burst from GRS 1747—312 which reached an es- 
timated maximum radius of almost 130 km. 

4U 1746—37, we used instead the estimated distances to the 
host globular clusters , Terzan 6 (9.5 kpc) and NGC 6441 
(11 kpc), respectively l|Kuulkers et al-lfeoOSt ). We found that, 
for all but one burst, / is strongly correlated with the maxi- 
mum radius (Fig. [3]). Excluding a single outl ier (an unusually 
inten se burst from GRS 1747—312; see also lin 't Zand et al.l 
2003), we find a Spearman's rank correlation of 0.742, sig- 
nificant at the 4.2<r level. 

We suggest that the preferential occurrence of high /- 
ratios in dippers as well as the correlation shown in Figure [3] 
can be accounted for either by a changing partial obscura- 
tion of the neutron star and its expanding photosphere by 
the disk during the radius expansion or by a varying reflec- 
tion of the burst emission off the far side of the disk. In 
the first case, which is illustrated in Figure [4] it is plausi- 
ble that even during non-dipping intervals, the neutron star 
in high-inclination systems is subject to a certain degree of 
obscuration and hence absorption by the disk material. The 
expanded photosphere present during the radius expansion 
episode increasingly protrudes above and below (relative to 
our line of sight) the disk, so that the emitted flux during 
radius expansion is less subject to absorption by the disk 
material, and the apparent flux is commensurably larger. 
When the photosphere touches down onto the neutron star 
again, the absorption increases once again, giving rise to a 
significantly lower touchdown flux. This requires a partic- 
ularly finely-tuned inclination, in which the viewing angle 
of the observer grazes the accretion disk at its maximum 
thickness. 

In the second interpretation, as the photosphere ex- 
pands, the far side of the disk intercepts and reflects an 




Figure 4. Schematic showing the suggested geometry giving rise 
to the low touchdown flux values (not to scale). Prior to burst igni- 
tion (top panel) the neutron star (dashed small circle) is partially 
obscured by the edge of the accretion disk. When a burst ignites 
and causes the photosphere to expand, the expansion out of the 
disk makes the photosphere more visible. As the photosphere col- 
lapses back towards the neutron star (lower panel) obscuration 
by the disk becomes increasingly important. 

increasing fraction of the emitted luminosity. This leads to 
an artificial increase of the observed flux at the peak of the 
radius-expansion episode, while the touchdown flux corre- 
sponds to the true Eddington limit. Albeit plausible, the 
large observed /-ratios require, in this interpretation, a large 
amount of reflection off the disk, which can by achieved only 
in the presence of significant warping. 

In either case, the ratio of the peak burst flux to the 
flux at touchdown is expected to increase with increasing 
radius of the photosphere, as indicated by the observational 
correlation shown in Figure 3. 

If the origin of the large /-ratios is related to interac- 
tions of the burst flux with the disk in the high-inclination 
sources, we might also expect a dependence on the orbital 
phase. Most of the dippers have only 1-3 radius-expansion 
bursts each; however, MXB 1659-298 and 4U 1916-053 
each have 12. We folded the burst times on the dip /eclipse 
eph emerides for these two sources l|Oosterbroek et alj|200ll . 
and lChou et al.ll200ll , respectively), and tested for variations 
in the /-ratio as a function of orbital phase. In neither case 
did we detect significant variation. In MXB 1659—258 the 
three largest /-values were all found in the phase bin 0.8-1.0 
(i.e. immediately prior to the eclipse); however, the variation 
was not significant compared to the scatter on the /-values 
in the other phase bins. 



4 DISCUSSION 

In this paper, we studied 246 photospheric radius expan- 
sion bursts from a large number of accreting neutron stars 
observed with RXTE to determine the systematic uncertain- 
ties that affect the measurement of their masses, radii, and 
distances using the touchdown (Eddington) fluxes. Compar- 
ing the peak fluxes to the touchdown fluxes of the bursts in 
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our sample, we found that for the majority of the neutron 
stars, this ratio is / < 1.6, which is expected on the grounds 
of the changing gravitational redshift during the radius ex- 
pansion episode (|Paczvnski fc Andersonlll986l ). However, for 
some neutron stars, this ratio can be as large as / ~ 7. We 
found that the only sources for which the ratio of these two 
fluxes is larger than / > 1.6 are dippers and Cyg X-2. 

For the dippers, we attribute the large /-ratios to geo- 
metric effects related to their high inclinations. Determining 
which of the two candidate mechanisms gives rise to the large 
ratios may be possible via more detailed X-ray spectroscopic 
analysis. For example, time-resolved spectral analysis of the 
radius-expansion bursts from dipping sources may provide 
evidence of a reflection component from the disk, although 
the necessity of deconvolving the changing contribution of 
this component from the underlying burst spectral evolu- 
tion likely presents a considerable challenge. Such additional 
analyses are beyond the scope of this paper. 

Cyg X-2 is the only source exhibiting bursts with large 
/-ratios (> 1.6) that does not also exhibit dips. This 
source, however, is distinct from the dip sources in sev- 
eral other respects. Cyg X-2 is the prototypical Z-source 
( Hasi nger fc van der Klisl 1 19891 ) and is thought to accrete 
at a near-Eddington rate, as inferred from obs ervations of 
its accretion flow in the UV (|Vrtilek et al.lll99Gu . This is in 
sharp contrast to the dippers, all of which accre te at sub- 
tantially sub-Eddington rates fe.g. lLiu et al .120071 ). The high 
accretion rate of this source may influence its behaviour in 
two different ways. 

First, we consider the possibility that the X-ray bursts 
observed from Cyg X-2 are of a diffe rent nature than t hose 
of the other bursting sources (see iTaam et al.1 1 19961 . for 
models of bursts in near-Eddington sources). Indeed, the 
bursts in the two high accretion rate sources, Cyg X-2 and 
GX 17+2, do not show, at first glance, the characteristic 
patterns of spectral evolution seen in the other bursters, 
making their classification as type I or typ e II bursts difficult 
jSztaino et al]|l986l ; rKuulkers et alll995l ). This discrepancy 
has been attributed, however, to artifacts introduced by the 
subtraction of the accretion flux, which is comparable to the 
burst emission in high-luminosity sources. 

Second, the large difference in the inferred accretion 
rates suggest that the accretion disk may be significantly 
thicker in Cyg X-2 than in the dip sources. As a result, the 
geometric effect discussed in §3.2 for the dippers may also 
be responsible for the large /-ratios found in Cyg X-2, if the 
latter is observed at a reasonably high inclination. 

Modeling the optical lightcurve of the bin ary constrains 
the inclination of Cyg X-2 to 50° < i < 75° (|Cowlev et al.1 
ll979l : IOrosz fc Kuulkerslll999h . Such large values of the incli- 
nation are also strongly suppo rted by the secular va riations 
in the X-ray flux of the source (|Kuulkers et al.|[l996T ) as well 
as by the super-orbital periodici ties found in its long-term 
lightcurve (|Wiinands et al.ll 1996h . The combination of a rel- 
atively high inclination with a geometrically thick disk per- 
haps makes Cyg X-2 similar to the dippers for the purposes 
of the observational effect reported here. 

The discrepancy between the maximum and touchdown 
fluxes for the high-inclination systems has implications for 
the measurement of neutron star masses and radii for burst- 
ing neut ron stars. F or example, the method recently pro- 
posed bv lOzell l|2006T i relies on the accurate measurement of 



the Eddington limit at touchdown. Because of this, it is, in 
principle, affected by systematic effects on the touchdown 
flux as the one presented here. However, the inferred val- 
ues of both the mass and the radius of the neutron star 
depend only on the ratio F coo \/Ft of the flux during the 
cooling tail of the burst to that during touchdown. Indepen- 
dent of the particular interpretation of the unusually large 
/-ratios observed in dippers, this ratio should remain un- 
changed because the bias depends on the size of the photo- 
sphere and hence affects in the same way both the cooling 
and the touchdown fluxes. The correlation presented in Fig- 
ure [3] further supports this conclusion. 

In contrast, we note that the distance D depends upon 
the ratio of the square root of the cooling flux to the Ed- 
dington flux, so that D is affected by the precise value of 
the Eddington flux. According to our first geometrical inter- 
pretation (where the touchdown flux is lower than the true 
Eddington flux due to geometric screening) the distance de- 
rived from the touchdown flux will be too large by a factor 
f 1/2 . For EXO 748-676, th is would indicate that the dis- 
tance derived by I6zell (|2006T l should be scaled by U (/) to 
yield a revised distance of 7.1 ± 1.2 kpc, using the mean ratio 
(/) = 1.7 ± 0.4 we calculate for this source. On the other 
hand, in the second geometric interpretation, the touchdown 
flux represents the true Eddington flux and the distance de- 
rived using this value is not affected. Thus, taking into ac- 
count both possible origins for the discrepancy between the 
peak and touchdown fluxes, the distance to EXO 0748—676 
is in the range 6-10 kpc. 

It is also worth considering other possible sources of er- 
ror on the mass, radius and distance derived via the touch- 
down flux and burst tail emission. In a previous study 
which attempted to determine the gravitational redshift 
from the variation of the burst flux du r ing p hotospheric 
radius-expansion episodes, iDamen et al.l fl990) concluded 
that systematic effects arising from variations in the (i) per- 
sistent emission, (ii) the shape of the burst spectrum, and 
(Hi) photospheric composition could not be neglected. The 
geometric effect described here illustrates that these three 
additional factors will only influence the derived mass and 
radius if they alter the ratio of the cooling flux to the Ed- 
dington flux, and the distance if they alter the ratio of the 
square root of the cooling flux to the Eddington flux. 

A possible change in the persistent emission during the 
burst may indeed affect either the touchdown flux or the 
cooling flux. Given that this effect will be most significant 
for sources with high persistent luminosities such as Cyg X- 
2 (as also shown by Damen et al. 1990), using low persistent 
luminosity sources to determine the neutron star mass and 
radius minimizes this uncertainty. 

As for the second effect, assuming the observed spectra 
are consistent with a blackbody (which is usually the case) 
the spectral shape will affect the derived fluxes only if the 
spectrum deviates significantly from a blackbody outside the 
instrumental bandpass. In the case of bursts observed with 
RXTE, the contribution of flux from above and below the 
bandpass is at most 10%, which likely represents a conser- 
vative limit on the degree of this effect. 

Finally, theoretical burst models as well as observations 
of 4U 1636-536 suggest that variations in composition may 
occur in the beginning of the burst, prior to the peak. In 
the touchdown and cooling phases, such variations are not 
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expected. Moreover, the color correction factors during this 
cooling stage, which enter the calculation of neutron star 
mass and radius, depend very weakly on metallicity (Ozel 
2006), thus limiting any additional bias. 
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